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Abstract 

In this second paper of the series on the coherence of wave packets, 
we study neutrinos in high energy experiments where neutrinos are 
produced by decays of pions or muons which are described using wave 
packets. The space time position where one neutrino is produced is 
not fixed to one value but is extended in macroscopic area. Hence the 
amplitude is defined by a superposition of the amplitudes of different 
neutrino's production time in the macroscopic region and depends on 
the absolute value of the neutrino mass. 

We analyze neutrino interference based on operator product ex- 
pansion near the light cone and find a new universal term in the 
time dependent neutrino probability. This new term has an origin 
in higher order quantum effect in a similar manner as axial anomaly. 
Roles of Lorentz invariance and the operator product expansion in the 
light-cone region are clarified and a possibility of measuring the abso- 
lute value of neutrino mass from neutrino interference experiments is 
pointed out. 
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1 Introduction 



One particle states in nature are described by wave packets in various situa- 
tions. In this series of papers [TJ [2j E] , we study physics and implications of 
wave packets. By using wave packets, space-time dependent amplitudes and 
probabilities that are impossible to obtain using standard scattering ampli- 
tudes are computed. The time dependent probability thus obtained is shown 
to give a valuable information on the absolute value of neutrino mass. 

Neutrinos are particles that are very light and interact with matters very 
weakly. Their masses were found to be finite from recent flavour oscillation 
experiments [U El El El E]- The neutrino oscillation experiments using 
neutrinos from the sun, accelerator, reactors, and atmosphere gave the values 
of differences of the mass squared. Their values are found to be [TO] 

Am 2 21 =m\-m\ = (0.8 ± 0.03) x 10~ 4 [eV 2 /c 4 ] (1) 
\Am 2 32 \ = \m 2 3 -m 2 2 \ = 0.19 to 0.3 x 10~ 2 [eV 2 /c 4 ], (2) 

with certain uncertainties where (i = 1 — 3) are mass values. The squared- 
mass differences are extremely small but the absolute values of masses are 
unknown. Tritium beta decays [UJ have been used for determining the ab- 
solute value but the existing upper bound for the effective electron neutrino 
squared-mass is of order 2 (eV/c 2 ) 2 and the mass is 0.2 — 2 eV/c 2 from cosmo- 
logical observations [T2]. Neutrino masses are far from other particle's masses 
and are important parameters of physics. Neutrino masses are not affected 
from standard model of electroweak gauge interactions and it is important 
to know precise values of neutrino masses. 

Neutrinos interact with matter by weak interactions and event rates are 
very low and neutrino detection is hard. However using its weak interactions 
with matters, neutrinos can be used as new observational means once the 
detection method is established [13] . Using neutrinos, several astronomical 
objects such as sun, moon, and other stars inside of which can not be observed 
directly by ordinary means such as lights, electrons and protons, would be 
studied in a future [2] . For these applications, it is necessary to know precise 
properties of neutrinos. We study wave and particle properties of high energy 
neutrinos. 

A double slit-like interference of a neutrino, which is totally different 
from the flavour oscillation, is a subject of the present work. Neutrinos are 
produced by weak decays of particles and propagate finite distance before it 
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is detected. The distance is not fixed but varies within certain range. So the 
wave at the detector is a superposition of the neutrino produced at different 
positions. We study the pion and muon decays and the neutrino detection 
probability in hadron collisions using wave packets. The scattering amplitude 
and probability are computed usually using plane waves but the finite time 
or finite distance behaviors are not known by the standard method and wave 
packets are suitable means to study these dependences. 

When the particles involved in a reaction are wave packets of finite spatial 
extensions and the particle measurement is made by wave packet, the time 
and space dependent transition probability can be observed. By using wave 
packets, the amplitude and probability of the finite time interval are also cal- 
culated. The probability obtained by our calculation has an oscillating term 
in time interval, T, in addition to the normal T-linear term. The normal 
term is calculable also in the standard S-matrix in the momentum represen- 
tation, but the anomalous T-oscillating term in the finite time is calculable 
only using wave packet. The wave packet formalism supplies the space-time 
dependent informations and we obtain the length dependent probability of 
finding neutrino. This dependence is generated since the position where the 
neutrinos are produced in the decay of hadrons is not fixed and extended to 
macroscopic area. The neutrino amplitude is a superposition of the ampli- 
tude of different time and position and shows the interference. 

We combine pion dynamics with the weak interaction of the neutrino and 
find the time dependent amplitude and probability. This amplitude shows 
neutrino interferences of anomalous behavior that is sensitive to the neutrino 
mass. The probability oscillates with time in a scale that is much shorter than 
that of flavour oscillation. The wave length of this oscillation is determined 
by the energy and mass of neutrino and the interference experiments could 
be useful for finding the absolute value of the neutrino mass. 

We have shown general features of wave packet scattering in [TJ and of 
particle coherence in the previous paper I [2] 0. An important feature of 
the wave packet of the relativistic particle is that the phase factor of the 
wave function is determined by the mass and the energy in a relativistic 
invariant manner. For the neutrino of mass m u and energy E u the phase 
factor is expressed by using the differences of two positions Ax = x — X 
and of two times At = t — T, as exp (i<j)), where the phase <ft is expressed as 
(f> = m v \/ c 2 At 2 — Ax 2 and also is written as <f> = x cAt, where (t, x) are 

1 The general arguments about the wave packet scattering are in [15l ITH1 ITT] . 
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time and space coordinates of the production point and (T, X) are those of 
the detection point. Consequently the interference due to this phase is also 
determined by the energy and mass. For the neutrino of very small mass, of 
order 1 eV or less, this interference is in the large scale. 

We investigate the physical problems that are connected with neutrino's 
wave packets and interferences in high energy regions. Particularly neutrinos 
from pion decay and muon decay are studied in this paper. Other low energy 
neutrino processes caused by solar neutrinos, reactor neutrinos, and others 
are studied in a next work. 

This paper is organized in the following manner. In section 2, wave packet 
sizes of decaying particles are estimated. In section 3, we study neutrino 
production amplitude and probability in hadron collisions and in section 4 
we study neutrinos from real pions decays and those of muon decay in section 
5. Summary and prospects are given in section 6. 

2 Wave packet sizes 

When decaying particles are not exact plane waves but are the wave packets 
of finite coherence lengths, the produced particles have also these properties 
of coherence. We estimate coherence lengths of proton first and those of pion 
and muon next following the methods of our previous works pQ [2J. 

2.1 Pion wave packets 

Pions are produced by the proton collisions with target nucleus. Hence the 
coherence property of pion is determined by the coherence property of proton 
and nucleus. The proton has a finite scattering probability with nucleus in 
the matter and has a finite coherence length and the target nucleus has a 
microscopic size of order 10~ 15 m and its position is fixed with the uncertainty 
of atomic distance in crystals. From these values coherence length of proton 
is determined and using the proton coherence length we estimate the pion 
coherence length. 

2.1.1 Proton mean free path 

As was shown in I [2], the coherence length of a particle in dense matter is 
determined by its mean free path. High energy particles in solid material 
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interact with atoms frequently and its mean free path is understood well. 
The mean free path is an average distance for one particle to move freely and 
maintains particle's coherence. A particle is expressed by one wave function 
in a finite distance defined by the mean free path. Beyond the mean free 
path, particles lose coherence and are expressed by a different wave function. 
Hence this particle state is different from the plane wave. The wave function 
which expresses this particle has a finite spatial size and finite momentum 
width. 

The mean free path of the charged particle is determined by its scatter- 
ing with atoms in matter by Coulomb interaction. The energy loss is also 
determined by the same cross section and is summarized well in particle 
data summary [10J. In order to estimate a mean free path of the proton of 
1 GeV/c, we use the proton energy loss rate for several metals such as Pb, 
Fe, and others as, 

- = 1 - 2 [MeVg^cm 2 ], (3) 
we have the mean free path of the 1 GeV/c proton 

proton = = n 9 , ln rj^-l 2 il = 50 " 100 M ■ ( 4 ) 

^Xj9 (1 — 2) x 10 [MeV g cm 2 g cm 6 \ 

At lower energy of 2 MeV/c, the energy loss rate is about 10 MeVg _1 cm 2 
and the mean free path is 

^proton = 10 [cm]. (5) 

Actually proton which is accelerated and reaches high energy at the end 
has a size of wave packet after acceleration, L a f ter , 

T -J ^after /p\ 

-^after — -^before x > [y ) 

^before 

where Lf, e f ore is the size before the acceleration and Vbefore and v a ft er are 
the velocities before and after the acceleration. The velocity is bounded by 
the light velocity c, and the ratio for reaching from the above momenta to 
10 GeV/c is from about 1.2 at 1 GeV/c to five at 0.2 GeV/c. Hence the 
proton of 10 GeV/c has the mean free path 



^proton 



40 - 100 [cm]. (7) 
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2.1.2 Pion mean free path 

Pions are produced by a collision of the proton with target nucleus. Coher- 
ence lengths of pions which are produced in high energy proton collisions with 
target nucleus are obtained using the above initial proton coherence size and 
target size. In relativistic energy region, particles have light velocity. Hence 
in the pion production, the coherence length of the pion 5xf is given from 
that of the proton 5xi as, 



5xi 5x f 
Vi v f ' 

Sxf — —5x{ Rj 5xi. 



(8) 



Consequently from Eq.([7j), the pion's coherence length of the momentum 
1 GeV/c or larger momentum is given by 

£ P ion ^ 40 - 100 [cm]. (9) 

We use these values of Eq. (JTj) and Eq. in latter sections. 

2.2 Muon wave packet 

The muon is produced from pion decay. By the decay of pion of finite coher- 
ence length, a finite coherence length of produced muon is determined. 

2.2.1 Decay of pion 

Coherence lengths of the muon is connected with that of the pion by the 
ratio of velocities, 

^•£pion '^•Z-muon (10) 



^pion ^muon 

and is expressed as 

r _ Vraxion r /-. -. \ 

"•^muon ^ L,J 'pion- V / 

^pion 

For the relativistic particles the velocities are light velocity and the velocity 
ratio is unity. 
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Since the initial pion has an extension of momentum Ap pion the final 
muon has also an extension of momentum Ap muon , 



APmuon = Ap pion + 0(A). (12) 



2.2.2 Muon coherence length 

Combining Eq. ([9]) and Eq. (TlTTl . the coherence length of muon is given by 

L muon « 40 - 100 [cm]. (13) 



2.3 Neutrino wave packet 

The size of wave packet for observed neutrino is determined by the object 
that neutrino interacts in detectors. Neutrinos interact with nucleus or with 
electrons in atoms. The nucleus have sizes of order 10 -15 m and the elec- 
tron's wave functions have sizes of order 10~ 10 m. We study the neutrinos 
described by the wave packets of these sizes in many particle processes. In 
this respect, the neutrino wave packet of the present work is different from 
some previous works of wave packets that are connected with flavour neu- 
trino oscillations jTH [191 1201 IZH [221 [231 121], where one particle properties 
of neutrino at production are studied. It is important to study the neutrino 
wave packet at the detector to study the interference. 
The muon neutrino interactions in detectors are 

v„ + e" e~ + (14) 

+ e~ -> \T + u e (15) 

z/ M + A + (A + 1) +X (16) 

v^ + A^v^ + A + X (17) 

The neutrino wave packet a u in processes (fl4j) and (fToT) is of order 10~ 10 m 
and the neutrino wave packet a u in processes (1161) and (j!7p is of order 10~ 15 m. 
In the following sections are discussed the neutrinos in short or intermediate 
baseline experiments. We will see that the neutrino production amplitudes 
are unchanged even though the smaller wave packet, 10~ 15 , is used. The 
reason why the result is unchanged even with such a small wave packet is 
that the neutrino is so light that its velocity v u is almost the light velocity. 
Consequently, the two space time positions of the neutrino are almost on the 
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light cone where the dominant contribution in the amplitude comes from, 
as it will be discussed in the next section. In fact the neutrino of energy 
1 GeV/c 2 and the mass 1 eV/c 2 has a velocity 

v/c = 1 - 2e (18) 

e= (^)2 = 5xl0 -19 j 

hence the neutrino propagates a distance 

I = Z (l - e )=l - 51, 5l = l x 5, (19) 

when light propagates the distance Iq. This difference of distance, 51 becomes 

51 = 5 x 1(T 17 [m]; l = 100 [m] (20) 
51 = 5 x 10~ 16 [m]; l = 1000 [m], (21) 

which are much smaller than the sizes of the above wave packets Eqs. (J9]) 
and (I13p . Hence, the neutrino amplitude at the nuclear target or the atom 
target should show interference. The geometry of the neutrino interference 
is shown in Fig. [1] 

The electron neutrino interactions in detectors are 



u e + e~ ->■ e~ + u e (22) 
v e + A -> e~ + (A + 1) + X (23) 
v e + A^e + A + X. (24) 

The neutrino wave packet a v in processes (122]) is of order 10~ 10 m and the 
neutrino wave packet a v in processes (123]) and (124]) is of order 10 -15 m . They 



are treated in the same way as the neutrino from the pion decay. 

Low energy neutrinos of order MeV from the sun or reactors must be 
treated separately and will be studied in a next paper. 

3 Neutrinos in hadronic collisions 

Applying the wave packet formalism, we obtain the space-time dependent 
neutrino amplitude. From this amplitude we study the transition probability 
at finite time interval. 
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Fig. 1: The geometry of the neutrino interference experiment. The neutrino 
is observed by the detector at T and produced at t\ or t 2 . 

3.1 Semileptonic weak Hamiltonian and decay ampli- 
tude 

Semileptonic decay of pion is described by the weak Hamiltonian 

H w = g J dxdfj,(f)(x)Jy_ A (x) = —igm^ J dx(j)(x) J§{x) (25) 

j v-a( x ) = Pi^'fO- ~ TsMs), M x ) = H x )0- ~ 1fs) u ( x )> ( 26 ) 

where (p(x), fi(x), and u(x) are the pion field, muon field, and neutrino field 
without any ambiguity. In the above equations, g is the coupling strength, 
ir(x) , Jy_ A (x), and Js(x) are the pion field, leptonic charged vector current, 
and leptonic pseudoscalar. 



3.2 Neutrino production amplitude in hadronic colli- 
sions 

Neutrino production amplitude from an initial state |aij) to an final state \fif) 
which includes the lepton pair is given by 

T = J tfxifiAHvWon). (27) 
9 



When the muon is chosen as a plane wave and other states are wave packets 
specified by central values of momenta, coordinates, and times, 



= IPproton; Nucleus) (28) 
= \p,k\v,k v ,X v ,T v \pn,X*,Tn,hadrons). (29) 

Expressions of particle states of the pion which ignores the life time and of 
the neutrino and muon are expressed by one particle's matrix elements 

(O|0(x)|p P i on ,X pion ,Tpi on ) (30) 

~ jV e~^^^^" T ^ 7r ^ TlT ^ e - iE (Ppion){t~T pion )+ip pion -('X-X pion ) 

(fx, fc; v, k v , X v , T v \p(x)^f 5 v{x)\0) (31) 



Xe i(B(^)t-p M .2) e i(s(p v )(t-T 1( )-p v .(2-^)) 

where the spinor's normalization is 

J2Hp,sMp,s)) = ^±^. (32) 

In the above equations, (130]) and (!3lil . 0"^ ,0"^ are sizes of the pion wave 
packet of the initial state and of the neutrino wave packet of the final state. 
To get the final expression of Eq. (130j) . the pion momentum p is integrated and 
is replaced by the central value p P i on - The pion wave packet was estimated in 
the previous section and the neutrino wave packet size is determined by the 
experimental apparatus. A neutrino interacts with a nucleus in the detector 
and the collision with nucleus occurs incoherently. Hence the wave packet 
size of the neutrino should be the nuclear size that is the unit of detector. 
To study neutrino interferences, we use the nuclear size for o v . 

The amplitude T for one pion to decay to the neutrino and muon is 
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written as, 



N' 



T 




(2tt)I \E(p,)J 



{Pf\<Kx)\<*i) 




In Eq. (133]) . (t, x) is the space-time coordinates where interaction takes place 
and the state fi'j is the hadron state in which one pion is taken away from the 
state (3f. From the dependence of the phase factor on the momenta of the 
muon, neutrino, and pion, by using integration by part it is shown that the 
average value of the energy and momentum is conserved. The integration of 
the momenta multiplied by the integrand of Eq. (1331) satisfies 



We use these relations later. 

Due to the coordinate dependence of the wave packets, the amplitude 
Eq. f )33|) depends on the space time coordinates. Furthermore, the integrand 
of Eq. f )33|) depends upon the space-time coordinate of the weak interaction 
where the neutrino is produced and is given as 



with a suitable normalization constant N". This amplitude depends upon 
the coordinates (t, x) explicitly and is not invariant under the translation. 
So this satisfies peculiar properties of translational non-invariant amplitude 
and the states of wide momentum region play the role. Even the infinite 
momentum states couple with <p(x) and appears in (3f and gives important 



(Ppion) = (Pf, + Pu) 
(Pjion) = ((P/l +Pvf)- 



(35) 

(36) 




(38) 



(37) 
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contribution to the probability at two different positions at finite times. This 
is quite different from the ordinary scattering amplitude where the infinite 
momentum state decouples from the final state due to the energy momentum 
conservation. We will study this point in detail later. 



3.3 Integration of neutrino momentum 

We compute the neutrino momentum integral of Eq. (13"3"|) in two methods. 
Gaussian integral is applied in the first one and stationary phase approxima- 
tion is applied in the second one. The former method is valid in small time 
interval and the latter one is valid in large time interval. In both methods, 
we obtain qualitatively same results. Especially the phase of neutrino wave 
function has a particular form that is proportional to the square of the mass 
and inversely proportional to the neutrino energy. 

3.3.1 Gaussian integral 

For not so large t — T u , the neutrino momentum p u integration of Eq. ( )33l) is 
made by Gaussian integral around the central momentum k v . The amplitude 
becomes then, 



igm^N I ^f(/3 / |0(x)|a i )e i(£{ ^ )t -^' x " ) M(p At )7 5 z/(4)e^ 

(39) 




where N is the normalization factor, v l u is the i-th component of the neutrino 
velocity, and <fi is the phase of neutrino wave function. They are given by 

h '" (f> = E(k v )(t-T v )-k v -(x-X u ). (41) 



V 



The phase factor of the neutrino wave function, 0, is rewritten by substi- 
tuting the central value x of neutrino's Gaussian function 

x = X v + v v (t-T v ). (42) 
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We have the phase 



E{k v ){t-T u )-k v -{x-X v ) 
E{k u ){t-T u )-k v -v u {t-T v 
E v \k v )-kl 



(43) 



E V {K 



which has a typical form of the relativistic particle. The phase becomes 
proportional to the neutrino mass squared and inversely proportional to the 
neutrino energy. 

3.3.2 stationary phase approximation 

Integration on the neutrino momentum p v in the amplitude at macroscopic 
time difference t — T u is made by stationary phase method. The stationary 
momentum is obtained from the stationary condition, 







dp 



~<P{Pv 



0. 



(44) 



The equation for the stationary phase Eq 
Pi 



E(p v ) 



- (0) 
Vv=Vv 



is written as, 

{t - t v ) - {x - xy = o. 



(45) 



The square of the above equation leads 



Pv 



E\p h 



Pv=Pu 



(0) 



1 + 



m 



Pv 



(t ~ %f 



Vv=Vv 



(0) 



[x-Xf 



(46) 
(47) 



Hence the phase factor at the stationary point is given by, 



Hpv% v=p ^ 



m[{t-T v y -{x-x v y 



1/2 



(48) 



m 



Ey{p v 



- - (0) 
Pv=Pv 



(t-T v 



This phase agrees to Eq. ( )43l) obtained by the Gaussian integral method. 
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Hereafter we use p v instead p} and we understand that the p^ in the 
following equations stands for p^ which is a function of time and space 
coordinate. The decay amplitude becomes 



When the space time coordinates (xq,x) are integrated, the delta function 
of the energy and momentum conservation appears. The scattering ampli- 
tude with this delta function has final states that have the same energy and 
momentum with the initial states. On the other hand, the space and time 
dependent amplitude T(t,x) is not invariant under the translation and has 
no delta function. So the energy and momentum of the final state is not nec- 
essary the same as the initial state. The states which do not satisfy the total 
energy and momentum conservation should be included to get consistent re- 
sults from the completeness. This scattering amplitude shows the space and 
time dependent behavior, which is a new information. So by interchanging 
the order of the integration, we are able to obtain the probability and other 
informations at the finite time interval. 

3.4 Decay probability and interference 

We find the probability of observing neutrino at certain finite distance. Due 
to the finite distance effect or finite time interval effect, the probability at 
the finite time interval is not invariant under the translation. So special care 
is needed due to the non-standard nature of the non-invariant probability. 
The integration over the momentum should be made with special care for 
the translational non-invariant probability. The infinite momentum states 
should be included from the completeness of the physical space unless en- 
ergy momentum conservation restrict. We will see that the state of infinite 
momentum actually contribute in a manner that is almost equivalent to vir- 
tual state. These pseudo-virtual states of the infinite momentum are actually 
important for the probability of finite time to get the consistent result. We 




igm^N / d^f(/3 / |0(x)K)e i ( E ®') t -^^iZ(p , , i )75i/(p I/ ) 



(49) 



(50) 
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compute the probability with two different methods which are applied in 
previous subsection. 



3.4.1 Gaussian integral 

Transition probability is a square of the above amplitude and is given by 



\T\ 



,2 2 



9 m a i I 

u n / 



47r\ 



N\ z I d A x 1 d A x 2 



x e i(E{p^-p^) x e -i(s(? M )t 2 -? M -3 a ) 



x e 2a 



(51) 



where S^Si, s 2 ) stands for the products of Dirac spinors and their complex 
conjugates, 



S 5 (sx, s 2 ) = (w(p At )(l - 7s)^(pl)) («(Pu)(l - ls)v{p v )y 



and its spin summation is 



S 5 = ^5 5 ( Sl , S2 ) 

-2(p M -Pv). 



(52) 



(53) 



1 



We use the relation Eq. ( 136|) and write S5 as 

1 



-(m pion -m M ). 



(54) 



3.4.2 Muon momentum integration 

When the muon in the final state is not observed, the muon momentum, p^, 
is integrated. Since integral on coordinates (t, x) is made later, the energy- 
momentum conservation does not hold for (t, x) dependent quantity. Hence 
the integration region of p M is whole momentum region from the completeness 
of the state. The muon of the infinite momentum is produced as a real state 
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but the time dependent amplitude should include this momentum region. 
Otherwise Lorentz invariance does not hold and meaningful result is not 
obtained. 

Let the following function be 

St = t 1 - t 2 , Sx = x 1 - x 2 , 

where time and space coordinates are fixed and so the four dimensional mo- 
mentum p M is not restricted from the energy momentum conservation. The 
integration region includes the infinite momentum where the velocity be- 
comes the light velocity. Consequently this function has a singular function 
near the light cone as 

AJSt,Sx) = -Ltf(A)e(tff) - i^=0(\) [/Vi (m.Vx) -ie(<ft)Ji (m^Vx) 

47T 87TVA l v > v 

+ i0(-X)-^= Kl (m.vQ) (56) 
A = (St) 2 - (Sx) 2 , (57) 

e(St) = W-«(-«) = f + ;f r ^n (58) 

1—1 for St < 

and gives the most important contribution to certain probability. We in- 
vestigate the effects of this most singular term later, which becomes most 
important in high energy limit and in loop amplitudes where the infinite mo- 
mentum states give finite contribution. Actually A M has several less singular 
terms and oscillating terms. These terms are not important and we ignore 
in the present work. 

The probability for not so large t — T v is computed from the amplitude 
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Eq. (I3"3"j) and is given by 



y muon 



|T| 2 (59) 



S1,S2 

3 3 

/47r\ 4 / 47r\ 4 

iV™=(— J - , 5t = h-t 2 , 8x = x x -x 2 . (60) 
3.4.3 Integration of the interaction position 

Next we integrate the coordinates x\ and x 2 . Integration on coordinates x\ 
and x 2 are made by the Gaussian integral on the neutrino wave packet and 
the A^St^x), 

Jo = J ^^e-^^-^^t' 1 -^)^-^^-^^^-^)^^^^^. 

(61) 

By changing the variables X\ and x 2 to X = Xl + X2 ; x = Xl ~ a:2 , the above 
integration is made easily, and we have 

^ifel*'-^ (62) 

Finally the total probability is written as 



y rfPmuon ^ |r| 



2 



(63) 



S1,S2 

2 



xe -^(^-X 7r +(v»-Vn)(t 1 -T v )+v 7T (TK~T v )) 2 e - 1 l_(x„-X w +(^- 1 / vr )(t 2 -T„)+^(T^-T„))" 
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From the pion coherence length obtained in the previous section, the pion 
Gaussian parts are regarded as 



, 2 



e -2^(^-^+(^~ t? -)( il - T -)+^( T -- T -))" ~ i (64) 

in a distance of our interest which is of order few 100 m. In a larger distance, 
this condition is not satisfied and the interference disappears then. This 
condition is that the neutrinos produced in the different decay area overlap 
each other. Other situations where this condition is not met, the interference 
pattern becomes different. 



3.4.4 Pion momentum integration 

It is convenient to classify the diagram for the fixed values of X\ and x?. of the 
Eq. (l5Tj) into two types. In the first one, the correlation function of two pions 
of coordinates X\ and X2 is composed of those of finite momenta of on mass 
shell and we write this as Fig. |2j In the second one, the pion correlation 
function is composed of those of infinite momentum and we write this as 
Fig. |3j Because the coordinates x\ and X2 are fixed, momentum and energy 
flows at these points can have arbitrary values in these amplitude, which is a 
peculiar property of position dependent amplitude. The infinite momentum 
flow from the X\ to the hadron part make some states in the final state (3 to 
have infinite momentum. Hence the final state of the space-time dependent 
amplitude, i.e., amplitude of fixed coordinates x\ and X2, should include the 
infinite momentum sates from the completeness of the particle states. This 
infinite momentum states contribute only to (t, x) dependent quantities and 
to the probability at the finite time. This states do not contribute to the 
standard scattering amplitude defined as the overlap between the states at 
t = ±oo. Hence this contribution is similar to higher order quantum effect 
or one-loop effect where virtual states of the infinite momentum give a finite 
contribution. We will see that this term gives the important contribution to 
the interference term at finite time. If the integration over the space time 
(t, x) is made first as in the ordinary scattering amplitude, the probability has 
no dependence on the space and time and these states of infinite momentum 
decouple from the integrated total amplitude at the infinite time. 

Pion in Fig. [2] is constrained by the initial states a and the final states 
/3 of the hadronic part and has the maximum momentum allowed from the 
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Fig. 2: Diagram of the normal term 
of transition probability for observing 
the neutrino produced at x\ or x 2 - All 
hadrons have finite energy and mo- 
mentum. 



Fig. 3: Diagram of the anomalous 
term of transition probability for ob- 
serving the neutrino produced at x\ or 
X2- One pion along one line has the 
infinite energy and momentum. 



energy and momentum conservation of the final state. The pion in Fig. [3j 
however, is not constrained and can have infinite momentum. Thus the 
two diagrams are different each others. However there is no way to separate 
directly two components using the observables. In stead of a direct separation 
of two components, the probability at the finite time is decomposed into two 
components, one is the T-linear term and another is T-oscillating term. Fig. [2] 
contributes to the the T-linear term and Fig. [3] contributes to the both terms. 
Consequently the T-oscillating term from Fig. [3] is uniquely computed from 
our formulation and we study this term. The pion momentum is taken in 
whole momentum region from the completeness of the states. The result 
gives a Lorentz invariant term of the variable A = (x\ — X2) 2 . 

The pion correlation function of space time coordinates is obtained by 
integrating the pion momentum and has two terms, 

A pion (5t,5x) = {P f \(f>(x 1 )\a i )(a i \<l>{x 2 )^ f ) (66) 

-^universal ( A) ~~\~ -?nomral (<^? <fe), 
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Fig. 4: The energy contour in the integration along the real axis is changed 
into that along the imaginary axis by Wick rotation. 

where the non-invariant term is given by 

F normal = P(p w )e-< E ^ Stl -^ sgl ) , (67) 

where the momentum in the above equation is that of the original pion 
momentum but not the momentum of the center of wave packet and is in 
the finite range allowed by the energy and momentum conservation. This 
correlation function is a function of the combination dx^P^-^ and is not 
invariant under Lorentz transformation of Sx only. The first one is manifestly 
invariant under Lorentz transformation of 5x and a function of A = (Sx) 2 . 

First term in the right-hand side of Eq. ( 166]) is due to Fig. [3] where the 
intermediate pion has the infinite momentum and the second term is due to 
Fig. [2] where the intermediate pion has finite momenta allowed from the en- 
ergy and momentum conservation. We argue further on the A pion (St, Sx) and 



A muon (St, Sx) from operator product expansions [25] in Appendix B Particu- 
larly the magnitude of .Funiversai(A) is estimated based on the Euclidean metric 
integration which is obtained by deforming integration path of the the in- 
termediate energy along the real axis into the one along the imaginary axis 
by Wick rotation as in Fig. HJ The integrand has a cut along the real axis 
due to infrared divergence and the integral thus defined is Lorentz invariant 
and would be valid as a regularized value. We estimated the integral in the 
Euclidean metric and found a finite value. We use this value in this paper. 

The non-invariant part comes from real pion decays where the pion mo- 
mentum has an upper bound. The invariant part comes from interference 
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terms where the pion momentum reaches infinity. We see that the amplitude 
from the non-invariant term oscillates with time as 

-^normal j /normal^ j (68) 

of a frequency determined by the pion energy as 

2 

77? 

= ^- (69) 

This frequency is a magnitude of a microscopic angular velocity that is much 
larger than 

The invariant term at A = 0, -F universa i(0) gives a dominant contribution 
at the light cone region. Hence, we substitute 

^pion(^) -^universal (0) ~\~ /normal^ j (70) 

in the following calculations. 

3.5 Probability 
3.5.1 Interference term 

Using the correlation function we have the final expression of the probability 
of observing neutrino at the time T or at the distance L = cT, 

dp muon y^\T\ 2 = N dt l( it 2 — — (1 + /ao rma ie^ (tl -* 2) ) (71) 

J rr„ Jo l*i - l v 

[-^universal (0) X C?(T, U>„) + /normal X gf(T, U T + 



S1,S2 



iVp ro b 



iV prob = (2v / 27r)V,/ M , f^g 2 ml(ml-ml), 
g(T,u) = t[ dX am{ ? X) + -{cos(o;T)-l}, (72) 



u v — -£T, w vr — -tt", — CI, 



where L is the length of decay region. Eq. (ITT]) depends on the neutrino 
wave packet size o v and the pion multiplicity but the number of neutrino 
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Fig. 5: The universal term of the neutrino probability at a finite distance is 
given by circles and those that subtracted the T-linear term shown by a solid 
line is given by triangles. The horizontal axis shows the distance in [m] and 
the probability is of arbitrary unit. The neutrino mass is m v = 1 eV/c 2 



events in the experiment is proportional to the initial energy, the neutrino 
reaction rate, the detector efficiency, and other parameters of the experiment 
in addition to Eq. flTTj) . We combine them to one constant iV ex p, then the 
event number of neutrinos at the distance L is proportional to 



Eq. ( 1721) has a slowly oscillating term and rapidly oscillating term. The 
latter term oscillates so rapidly that in the ordinary experiments only its 
average is seen. Hence, this term is regarded as a linear function of T, 



exp 



(-^universal (0) X g(T, U v ) + /normal X g(T, U n + L) u )). 



(73) 



E 



V 




(74) 
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and we have the probability 



/ dp muon ^ \ T \ 2 = — [-^universal (0) X g(T, U u ) + /normal X T] . (75) 

The probability of observing the neutrino at the time T is plotted in Fig. |5]for 
m u = leV/c 2 . From this figure we see that it consists of the T-linear term and 
the small and slow oscillation term on top of the T-linear term. The T-linear 
term comes from the invariant term and also from the non-invariant term. 
So the relative magnitude of oscillation term is not definite but from the 



calculation presented in |Appendix B[ we expect that the magnitude of the 



universal term is about the same as the the non-universal term. Because the 
probability at the time T has a linear term and the universal slow oscillation 
term, the universal term is extracted by subtracting the T-linear term from 
the total probability. The slowly oscillating term thus obtained is given in 
Fig. |6j This shows the slow oscillation of showing small neutrino mass. Thus 
we find the final form of the probability which has the oscillating term with 
the frequency determined by the neutrino's mass and the energy and T-linear 
term which comes from the Lorentz non- invariant term in Eq. (I66p . 

The slowly oscillating term is determined by the Lorentz invariant term 
in Eq. fl66|) . A typical oscillation length Lq is 

L ° H = " 2 °°^eW' (76) 

This length is quite interesting value for the observational point. We use 
m v = 1 eV/c 2 throughout this paper. 

The oscillating probability depends upon the energy E u and time T and 
its signature is observed also in the energy dependence of the probability 
at a fixed T. This total probability multiplied by the neutrino energy is 
plotted in Fig. [7] and the the probability obtained by subtracting the smooth 
background from the total probability is given in Fig. EJ The slow oscillation 
due to the universal term is seen clearly. 

3.5.2 muon in pion decays 

When the muon is observed in the same processes, the anomalous oscillation 

2 

is determined by the muon mass and energy as and the muon mass is 
larger than the neutrino mass by 10 8 . Hence the oscillation length is smaller 
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Fig. 6: After the T-linear term is subtracted, the probability shows clear 
oscillation. The horizontal axis shows the distance in [m] and the probability 
is of arbitrary unit. The neutrino mass is m v = 1 eV/c 2 

than that of the neutrino by 10 16 . For the muon of energy one GeV/c, the 
oscillation length is order 10~ 12 m. This value is too small to observe in 
experiments. 

Since it is hard to see anomalous oscillation of this length, it is meaningful 
to study an average probability. The total muon flux includes the effect of 
the anomalous term and is written as 

F ln = C(n w ) (78) 
F la = C'(n ir ), (79) 

where F\ n is due to naive decay of pions without interference effect and so 
is proportional to pion number, whereas i*\ a is due to the interference term 
that gives the anomalous oscillation of the neutrino. Because we ignored the 
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Fig. 7: The energy dependence of the neutrino probability is plotted. The 
probability at T is given by circles and those that subtracted the T-linear 
term shown by the solid line is given by triangles. The horizontal axis shows 
the neutrino energy in [GeV] and the vertical axis shows the neutrino proba- 
bility in the arbitrary unit. The neutrino mass is m„ = 1 eV/c 2 . The distance 
is 200 m 



non-singular term in Eq. (tool) , it is hard to know the magnitude of rapid 
oscillation term and the total flux. We expect that the Fi a has a same 



magnitude as F ln from the calculation of |Appendix B 



3.5.3 Stationary phase approximation 

So far the calculation is based on the Gaussian integration on the neutrino 
momentum. The result is the same in another calculational method of using 
the stationary phase approximation. 
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Fig. 8: The energy dependence of the oscillating neutrino probability that is 
subtracted the background term from the total probability is plotted. The 
horizontal axis shows the neutrino energy in [GeV] and the vertical axis shows 
the probability in the arbitrary unit. The neutrino mass is m v = 1 eV/c 2 . The 
distance is 200 m 



The transition probability is a square of the Eq. [49] and is given by 



\T\ 



^] 2 \N\ 2 I d A x l d i x 2 



5*5(51,52) 



m,. 



E 1 E 2 

V V 



,2 



x e ^(* 1 -^) e -<^(t 3 -r„) e -^((^-fc„) a +(^-fc„) a ) 

x e -i(E(p p ion)(t 1 -T 7r )-p pion ix 1 -jC lr )) x e i(E(p pion )(1?-T*)-p p i on -(&-jt*)) 



x e 



x e 2ct t 



(80) 



where Ss(si, sz) stands for the products of Dirac spinors and their complex 
conjugates defined in Eq. fl52|) and its spin summation of Eq. fl54|) . The prod- 
ucts of 5*5 with pion momentum is computed similarly and the probability 
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-'muon 

S1,S2 



for the unseen muon in which muon momentum is integrated is 

Em 2 (si) 

S1,S2 

= \N^\ 2 J d 4 x,d A x 2 |r^ )Qi | 2 (-^) 2 e -^-^ + ^-^%-< E ^^-^ 5i ) 

The final expression has a simple form and is the same as the Gaussian 
integral methods. 

4 Neutrino in real pion decay 

So far the initial state are the proton and nucleus and the pion is unidentified 
and is equivalent to intermediate state. 

We study the case where the pion is identified and is in the initial state 
of the scattering. Let the momentum of the pion is g P i on , then the decay 
amplitude of the pion is given as 

T = gm a J d 4 x (0\<j>(x)\q pion ){p^iy(x)\J 5 (x)\0) (82) 

= gm^J d 4 x e- i(q ^- pvyx (p^\fi(x)\0)(l - -f 5 )u(p u ) 

After the integration of the momentum, the total probability becomes 

dp m uon\T\ 2 (83) 

g 2 ml J d 4 Xl d 4 x 2 e-^^-P^-^A^X! - x 2 )S 5 



In the above result, the oscillation is so rapid that can not be observed and 
the time average is found as 



'muon | 



T| 2 (84) 
g 2 m\{ml - mj) J d A x 1 d 4 x 2 e - i(q ^~ PvHxi ' X2) A^ - x 2 ) 
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Fig. 9: The neutrino probability of the anomalous term is shown by circles 
and those that subtracted the T-linear term (solid line) is shown by tri- 
angles. The horizontal axis shows the distance in [m] and the magnitude 
is of arbitrary unit. The small oscillation is seen. The neutrino mass is m u = 
leV/c 2 . 



The above neutrino probability has the same property as the normal term 
Eq. fl68|) and shows no long distance interference. 

In higher order effect, the momentum of the intermediate state includes 
infinitely large value, so the light cone behavior is modified. Consequently 
it is expected that the long distance interference is generated. The lowest 
higher order correction is given by the pion exchange term. In this case the 
external pion momentum is included in the momentum of the intermediate 
pion and its effect remains in the final result, 



\T\ 



15) 



g g W7r Id x x d x 2 A^(xi - x 2 )A 7r (p pion , x x - x 2 )S 5 e ip "' 



where the correlation function A^Qopion, x x — x 2 ) depends on the external pion 
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Fig. 10: The length dependence of the oscillating neutrino probability that 
is obtained subtracting the T-linear term from the total probability. The hor- 
izontal axis shows the distance in [m] and the magnitude is of arbitrary 
unit. The neutrino mass is m u = leV/c 2 . 



momentum and is not a function of A. A 7r (p P i 0n , x\ — x 2 ) is computed in the 



Appendix B and becomes 



(Ppion; 

1 



%2 



(86) 



1 1 



PpkmOi - x 2 ) 64?r 2 m 2 



The probability of observing the neutrino in the real pion decay is given 
in Fig. (£)]), which is composed of the T linear term and the oscillating term. 
Because it is hard to know the relative size of the one loop amplitude with 
the tree amplitude, the relative magnitude of the oscillation term is arbitrary. 
It is known that the non-universal terms from the tree level and the one-loop 
level behaves like the T-linear term and the universal term oscillates. Hence 
the oscillating term is obtained by subtracting the T-linear term from the 
total probability. The result is plotted in Fig. ( TTUj) as a slow oscillation. 
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5 Neutrinos from muon decay 



5.1 Leptonic weak Hamiltonian and three body decay 
amplitude 

Muon decays to an electron and two neutrinos. One is an electron neutrino 
and the other is an muon neutrino. We study the flux of each neutrino at a 
certain distance L in high energy region. 

Leptonic decay of muon is described by the leptonic weak Hamiltonian 



without any ambiguity. In the above equations, Gf is the Fermi coupling 
constant, ji{x) is the muon field, e(x) is the electron field, v^x) is the muon 
neutrino field, and u e (x) is the electron neutrino field. Jy_ A (x) is the leptonic 
charged current. 

5.2 Neutrinos from muon three body decay 

Muon decay is treated in a similar manner as the previous pion decay. Here 
we assume that the muon correlation function has two components. First 
one is a function of A = 5x 2 and is invariant under Lorentz transformation of 
the coordinates 8x. Second one is a function of the product Sx^P mnon ^ and 
has the same property as the normal term of the pion decay. 

So the decay amplitude and decay probability are computed in Gaussian 
integral method. We integrate unseen particle's momenta and average over 
the initial muon momentum. 

Leading term of expectation value of the weak leptonic currents at the 
light cone is 



Because the electron and one unseen neutrino are not observed and their 
momenta are integrated in the infinite regions and the above correlation ([59]) 
is more singular than the previous one particle case. 




(87) 



Jv-a( x ) = H(%)l a ( l - + e(x)7 a (l - j^)v e {x) 



(88) 



(OlJ^OJ^Or^lO) 

= 8(-d 2 W2 + dMie^xl - xl))5'{\). 



(89) 
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The normal term gives the rapid oscillation term and its time average 
is proportional to the time T. The one loop term is generated by QED 



correction. As is discussed in the |Appendix B[ the fully invariant term does 



not exist in QED one loop correction but the non-oscillating power correction 
term which becomes roughly a times the magnitude of the normal term, 

/universal Q! ~, r X I /"normal I (90) 

do exist. This term gives an oscillation of the frequency determined by the 
neutrino mass and energy. 

Finally we have the similar expression for the decay probability as the 
previous pion decay, 

J2\ T \ 2 = M9i(T,"v) + 92(T,u u )} + f,„ h (91) 



^ S1,S2 

9i(T,co) = T£dX^^l, g 2 (T,u) = j\x 



sin(o;X) frr N f ^ v sm(uX) 

~X ' 

Although the oscillating term that is due to the invariant term is much smaller 
than the ordinary term, the observation of this term may be possible. 

In muon decay one muon neutrino and one electron neutrino are produced. 
They are linear combination of three mass eigenstates and a unitary MNS 
matrix combines flavour eigenstates with mass eigenstates. 

5.2.1 Electron from muon decay 

If the electron is measured and two neutrinos are unseen we have the prob- 
ability of the electron at the distance L, 

J dp u dp u i ^ \ T \ 2 = N e {gi(T, LO e ) + g 2 {T, U e )} + / e ,n rmal, (92) 



,2 



Tfl 

= , (93) 

where the oscillation length is given by 

Li = (94) 
mi 
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The oscillation length of the electron is given by 

L e = 10" 8 [m] x (£ e /GeV) 



(95) 



and is too small if the energy is a few GeV. In the ultra-relativistic energy 
region, this oscillation may become observable. 

6 Summary and implications 

In this paper, we showed that one particle states of decaying particles that 
produce neutrinos are described using wave packets of finite coherence lengths 
and studied its implications to the neutrino interferences. 

The wave packet size was determined either from particle production 
processes and detection processes. In the former, a finite mean free path in 
matter is the origin of the wave packet. The finite mean free path makes one 
particle to have a finite spatial extension and a finite momentum uncertainty. 
The state of a finite mean free path is a non-stationary state and is varied 
with time and space. In the latter, a finite size of the unit of detector is the 
origin of the wave packets. The wave packet sizes of the proton, pion, muon, 
and the neutrino were estimated and were used in analyzing high energy 
neutrino reactions. 

Since the overall phase of wave packet during propagation is determined 
by the time component that is proportional to the energy and the space com- 
ponent that is proportional to the momentum and the space position and the 
time position are connected each others, both effects are taken into account 
simultaneously. Due to the relativistic invariance in the energy and momen- 
tum and in the time and space position, both terms in the total phase are 
almost cancelled and the total phase becomes small number that is propor- 
tional to the mass squared and inversely proportional to the energy. Conse- 
quently when the neutrino is described by the wave packet, the space-time 
dependent probability of the neutrino is found and the above overall phase 
or its difference becomes observable. We showed that the time dependent 
interference of the neutrino in the processes of the decays of pion or muon 
reveals this phase. 

The time dependent probability of observing the neutrino at finite dis- 
tance was calculated for high energy collisions and the anomalous oscillating 
term was found. This term has the origin in the higher order quantum effects 
where the infinite momentum virtual states play the important role. The new 
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universal term is manifestly invariant under the Lorentz transformation of 
the coordinates and gives the most important contribution in the operator 
product near the light cone region. Because the neutrino's velocity is almost 
the light velocity, the time dependent probabily of finding the neutrino is 
determined by this universal term. The probabily of finding the neutrino at 
finite medium time is oscillating with the slow angular velocity in Eq. (T75|) . 
Since the angular velocity is determined by the neutrino mass and energy, 
the absolute value of the neutrino mass would be found from the neutrino 
interference oscillations. 

Due to the relativistic invariance, the correlation function and others 
become functions of the Lorentz invariant combination A. The space-time 
points that satisfies A = are on the light-cone surface and infinite number 
of points are on the surface. This is a feature of a relativistic invariant system 
and is a reason why the interference of the present work occurs. For a non- 
relativistic system, in a stationary state of the same calculation of the space 
coordinates is made by, 



and the only one point 8x = satisfies the condition and the probability get 
a contribution from only the point 8x = 0. The rotational invariant three 
dimensional space is compact but the Lorentz invariant four dimensional 
space is non-compact. This difference is important for the reason why the 
relativistic system has a peculiar property of the interference. 

It is worthwhile to clarify the difference of the space-time dependent prob- 
ability of the present work with the normal scattering probability defined at 
t = ±oo here. The normal scattering amplitude is defined from the overlap 
between the in-state at t = — oo and out-state at t — oo, and the space and 
time coordinates are integrated from — oo to oo and the energy and momen- 
tum of the final state is the same as that of the initial state. Hence the 
momentum of the muon or the pion in the final state of the ordinary scat- 
tering experiments are bounded due to the energy momentum conservation. 
So the infinite momentum is not included in the muon or pion of the final 
state. However the amplitude and probability at the finite time and their 
behaviors at the finite time are not computable in the ordinary S-matrix. 

In our method it is possible to compute the amplitude and probability 
at the finite time and space. The energy and momentum conservation does 
not hold for these quantities and the infinite momentum state of the muon 




(96) 
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and pion are included. These states of the infinite momentum give the finite 
contribution to the time dependent probability but do not contribute to 
the cross section measured at infinite distance. The important informations 
are obtained from the wave packet formalism that are not calculable in the 
standard scattering amplitude. Hence our calculation does not contradict 
with the ordinary calculation of the S-matrix in momentum representation 
but has the advantage of giving new informations. 

In our calculation, Lorentz invariance is one important ingredient. 

The characteristic small phase of the relativistic wave packet shows macro- 
scopic interference of the neutrino. Although this result should be applied 
in high energy region, it would be interesting to see if this effect is found 
in ground experiments and others. Depending on the mass value, the phe- 
nomenon we have discussed in this paper may be relevant to short base line 
experiments, long base line experiments, and atmospheric neutrino experi- 
ments and others. 

The oscillation phenomenon of the present work is sensitive to small mass, 
hence the same mechanism would work if there exists a very light particle. 
A possible candidate of light particle is axion. Axion might show a peculiar 
oscillation if it exist. 

In this paper we ignored the effects of the pion life time and the pion mean 
free path in studying the higher order quantum effects. We will study these 
problems and other large scale physical phenomena of low energy neutrinos 
in subsequent papers. 
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Appendix A Formula for neutrinos of three 

flavour 

There are three eigenstates of neutrinos and their mass difference squared are 
known and mixing parameters are also known partly from flavour oscillations. 
The wave packets have been studied in flavour oscillations in [IS] [HE [2H1 ED 
|2"2"| [221 [21] • We studied neutrino spatial interference in this paper, which 
is unrelated directly with these flavour oscillations. A unified treatment of 
neutrino phenomena is possible and will be presented in a next paper. 

For the wave packets to overlap at the detector, two components of mass 
eigenstates should arrive to the detector same time and should have the 
same energy within wave packet uncertainties. Hereafter we assume that 
these conditions are satisfied and study the amplitude for three neutrinos. 
General finite time amplitude for each flavour combination is expressed using 
amplitudes of three mass eigenstates, T(i, ti,Xi) as 

T a ,i3= dtidxi U a jT(i, ti, x^Ujg, (97) 
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where C/j Q is a unitary matrix which connects flavour eigenstate a to the mass 
eigenstate i. Time dependent neutrino production amplitude T(i,ti,xi) was 
given before and is substituted into the finite time probability. We have, 
after the coordinates and other variables are integrated, 

l-T 

\T a , P \ 2 = u ail {u ah y dhdhTih^Tii^hyu^iu^y, (98) 

JTo 

1 2 
m ■ m ■ 1 

h)T{i 2 , t 2 )* = Ne'^ 1 -^-^^-^ , (99) 

^i — h 

where N is a constant. 

Appendix B Operator product expansion and 

new universal term 

B-I Pion decay 

Neutrino detection amplitude and probability we have discussed is under- 
stood from operator product expansions at the light cone region (xi — X2) 2 ~ 
0[25]. 

In Eq.f l5ip . the total probability is given by the integral of the space time 
coordinates of the weak Hamiltonian and is invariant under the translation 
of space and time. Due to the translational invariance, the energy and mo- 
mentum of the final state are the same as those of the initial state. Now we 
interchange the order of the summation of the coordinates and final states 
and obtain the time dependent probability. The time dependent probability 
which is obtained by summing the final state first is given in Eq.( l7"TT) . This 
probability has two components, invariant term and non-invariant term un- 
der the translation in time. The former is the T-linear term and the latter 
is the T-oscillating term. 

The translational invariant term gets the contribution from the final states 
that has the same energy as the that of the initial state. But the energy of the 
final states that contributes to the latter is not necessary be the same as that 
of the initial state. Hence the states of the infinite momentum could appear 
in the intermediate state for the T-oscillating term and give the finite contri- 
bution to the probability at the finite time although they do not contribute 
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to the probability at the infinite time. The states of infinite momentum con- 
tributes only to T-oscillating term and we estimate its magnitude based on 
the operator product expansion at the light cone region. 

The functions A muon (<5t, Sx) and A pion (5t, Sx) at the light cone region, 

A = (St) 2 - (Sx) 2 = (100) 

are the expectation values of the products of the muon field and the expec- 
tation value of the pion field 

(0|/i(xi)/x(^)|0), (101) 

(«in, Pont 1 4>n{xi) (pn(x 2 )\ ttin, /?out), (102) 

where a- in includes the proton and the target state and (3 out includes many 
pions and other particles in the final state. 

The muon correlation function was studied in Eq. (|56|) . 



B-I.l the normal term 

Hereafter we study the pion correlation function. The pion correlation func- 
tion from Fig. Ogives the neutrino production probability that is determined 
by incoherent decay of produced pions and interference among the pions are 

.m 2 

negligible. The time dependence is given by e*"^* 1- * 2 ^ and the angular veloc- 
ity is too large to observe the oscillation. The correlation function becomes 

A pion (St,Sx) = J dp P(p)e lI £^- t2) (103) 

= (\N n \)e l ^ {tl - t2 \ (104) 

and in the ordinary pion energy of about 10 GeV/c 2 , the angular velocity of 
this term corresponds to the oscillation length 

L = ^£ (105) 

which is too short for observation. The average of this term contributes to 
the T-linear probability. 
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Fig. 11: The Feynman diagram of the 
anomalous term of the neutrino prob- 
ability using quark fields is given. 

B-I.2 a new universal term 



Fig. 12: In the anomalous term of 
the neutrino probability, the hard pion 
couples either with the pion or the pro- 
ton in the hadronic parts. 



The energy of the final state (3 in the T-oscillating term is not necessary the 
same as that of the initial state and infinite momentum states can couple 
and gives finite contribution to the pion correlation function. We estimate 
the effects of the infinite momentum hereafter. 

We study the pion correlation function that includes infinite momentum 
states in (3. They are almost equivalent to Fig. [2] except that pions have 
infinite momenta and we write the infinite momentum pion explicitly as 
Fig. |3J We will see that this gives a manifestly invariant term F(X) of Eq.f l66p . 

This correlation function is calculated with quark fields as in the Feynman 
diagram of Fig. [11] or with pion fields as in Fig. [3j In the former, QCD is 
applied. The amplitudes in QCD has severe infrared divergence in time-like 
region and quark propagator has no simple pole and has a cut at real energy 
axis as in QED discussed next. We replace the integration on the energy in 
Minkowski metric to the Euclidean four momentum integration. This Wick 
rotation is allowed if the amplitude is analytic except cut along real axis. 
We assume that this holds and compute the integral in Euclidean metric. In 
Euclidean metric there are no infrared divergence and quark dynamics are 
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effectively described by meson dynamics. So we compute the pion correlation 
functions using pion propagators as in Eq. tppj) . 
Now we estimate the amplitude. 

F((x 1 -x 2 )) (106) 

d'pe'^-^ - \ , - L- — 5{(p + Pf y-m 2 )\Tl- I 2 



d 4 p e ~ ip < xl - X2) - 



1 1 

n 2 — iep 2 — rn- 
1 1 I 



"5 Q __ 2 o" i • VV.F ' Fj J 111 ; Khadron 



x f \\r l I 2 

V /„., ... \o _<? • / \o ™9 1 /I hadron I ' 



^2 _ m 2 _ ^ £ ^2 _ m 2 _|_ ^ 2i7T 

1 

(p + pf) 2 — m 2 — it (p + pf) 2 — m 2 + ie' 

where pf = p a — pp. The momentum of the pion which connect the coor- 
dinates X\ and x 2 must be taken from —00 to 00. The integral is reduced 
to 

F{{x 1 -x 2 )) (107) 
-IT 1 I 2 1 l — + G 

— 1 1 hadron I ' , _ . „ „ "I" ^ 



Pf(xi - x 2 ) 647T 2 m 2 \ {p f (x 1 - x 2 )Y 

for pf 7^ in the Euclidean metric space calculation. The infinite momen- 
tum state of the pion at the coordinate x\ enters the hadronic part 7h adron 
and interacts with the one of pions or proton and the pion of the infinite 
momentum goes out because other particles in the hadronic part have finite 
momentum as Fig. [T2l 

For pf = 0, the invariant term is computed as 

F{{x x -x 2 f) (108) 

2 f d *PE frg (»!-»„) ( 1 ^ | T 1 12 

5 y (2tt)^ vpI+^v hadron ' 

The integral in the right-hand side 

is invariant under the Lorentz transformation of X\ — x 2 and is a function of 
A. Hence we have 

(0) ~ 7 \pi + m 2 ) "cm^ 2 ^' (110) 
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and 



F(0) (111) 

= [dhp E ( i y lr i |2 

y / /o„\4 \ „2 i ,™2 7 I hadronl 



(2vr) 4 VpI + ^ 2 



1 



IT 1 



|2 



647r 2 m 2 ' hadronl ' 

Actually because two pion couples with the proton or pions and the total 
number of the pions N n is larger than the number of the proton. Hence the 
pion correlation function becomes 

F{0) = 64^ ^l^ronl , (H2) 

where g^N is the pion Nucleon coupling strength and g nw is the pion pion 
coupling strength, 



B-II Real pion decay 

When the real pion of the momentum p P i on is an initial state, one loop cor- 
rection to the pion correlation function is given by the following integral, 

A n (p pion ,x 1 - x 2 ) (113) 

= / d 4 qe-^-^- \— ] —-5((p + q pion ) 2 ) 

J p z — m A — iep z — m z + te 

= [ d 4 q c ~ iq ' {xi ' X2) 1 - — 

J ' p 2 — m 2 — ie p 2 — m 2 + ie 2in 

1 1 

X (7 ^ 9 -~1 ^ 9 -)■ (114) 

[P + q P ion>- m 2 - ie {p + q P ion> 1 - m 2 + ie 

When the higher order correction is added and the infrared divergence is 
avoided, the propagator becomes not to have a simple pole but a cut along 
the real axis. Although it is a difficult problem to find out the dynamics of 
infrared divergence, we simply make Wick rotation and compute the Lorentz 
invariant term. 

This integration is made in Euclidean metric space. We have then, 

^tt (Ppion; 

x 1 -x 2 ) (115) 

1 1 1 + o< 1 



PpionOn -x 2 ) QAn 2 m 2 \ {p pion ( Xl -x 2 )} 



2 
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Fig. 13: The diagram of the normal Fig. 14: The diagram of the lowest or- 
term in the neutrino probability in the der correction in the neutrino proba- 
muon decay. bility in the muon decay. 



which has a power term. 



B-III III muon decay 

For the muon decays, higher order corrections are generated from QED and 
the hard photon exchange term of Fig. [TD would give the invariant term. 
This diagram has the infra-red divergence that is avoided by redefining one 
particle charged state in the initial and final state in such way that is dressed 
by soft photon. Charged field becomes not to have simple pole but cut. So 
integration of the large momentum in Fig. [TU becomes equivalent to that of 
the Euclidean metric. We estimate the Lorentz invariant amplitude from the 
integration in the Euclidean metric integration and we have then 

(xi - x 2 ) (116) 

= 2 f d4 PE ip E (x 1 ~x 2 )( 1 \2 1 | T 1 |2 

J (2vr)4 >| + m^ (p E + qf 

e 2 1 1 ! 2 1 
~ 64 7 r 2 g(x 1 -x 2 )mj' muonl + \q{x 1 -x 2 )) 2)) 

where q is the momentum of muon. Thus higher order correction has a power 
correction term that is inversely proportional to the variable q(x\ — x 2 ) but 

2 

not a function of A. This term gives also the oscillation of the frequency r ^ L . 
However because of the factor -<— ^ — r this oscillation appears in the time 
derivative of the probability. 
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